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INTRODUCTION 

With  the  advent  of  nanosecond  lasers,  new  opportunities 
have  arisen  as  well  as  problems  of  scientific  and  military  Interest. 
One  problem  Is  the  measurement  of  the  polarization  state  of  nano- 
second laser  pulses.  We  solved  this  problem  by  developing  a novel 
four-channel  polarlmeter  which  simultaneously  measures  the  four 
Stolces  polarization  parameters  of  an  optical  pulse.  This  design  Is 
enhanced  by  proper  selection  of  the  polarization  state  of  the  trans- 
mitted beam,  which  makes  It  possible  to  determine  the  ABCD  polariza- 
tion matrix  of  an  optical  system,  using  a single  optical  pulse. 

The  polarlmeter  Is  automated,  has  no  moving  parts,  and  Is 
controlled  by  a Hewlett  Packard  (HP)  9825A  desk  calculator  and  an 
HP-IB  Interface  Bus  System  with  an  HP  6940B  multiprogrammer. 

Two  Important  factors  were  critical  to  the  adoption  of  this 
polarlmeter  design.  The  first  was  the  realization  that  If  an  optical 
system  was  Illuminated  with  a nanosecond  laser  pulse,  the  traditional 
methods  for  measuring  the  polarization  state,  or  equivalently  the 
Stokes  polarization  parameters,  are  completely  Inadequate.  In  order 
to  measure  the  Stokes  parameters,  four  distinct  settings  of  a wave- 
plate/polarlzer  combination  are  needed.  It  Is  obviously  Impossible 
to  do  this  In  such  a short  time  frame  when  dealing  with  nanosecond 
pulses.  Since  four  distinct  settings  must  be  made,  a polarlmeter 
must  contain  four  separate  channels  In  which  the  wave  plate/polarlzer 
Is  preset  and  permanently  fixed.  With  this  design  It  Is  then  pos- 
sible to  determine  the  Stokes  polarization  parameters  simultaneously. 
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We  wish  to  emphasize,  however,  that  this  arrangement  allows  us  to 
measure  the  Stokes  parameters  of  an  optical  beam  Irrespective  of  any 
time  duration  of  the  signal. 

The  second  important  factor  In  the  design  of  the  polarl- 
meter  was  the  fact  that  most  optical  systems  of  military  Interest 
contain  only  polarizing  and  phase  shifting  elements,  as  we  shall  show 
in  the  following  sections.  This  being  the  case.  It  Is  possible  to 
show  that  the  4x4,  sixteen-element,  Mueller-Stokes  polarization 
matrix  has  only  four  unknown  matrix  elements.  A,  B,  C,  and  D. 
This  specialized  matrix  form  is  called  the  ABCD  polarization  matrix. 
Further  analysis  shows  that  If  an  optical  component  or  system  Is 
Illuminated  with  either  circularly  polarized  or  linearly  polarized 
light,  each  of  the  detected  Stokes  parameters  Is  equal  to  one  of  the 
A,  B,  C,  or  D matrix  elements. 

The  design  and  operation  of  the  polarimeter  are  directly 
related  to  the  polarizing  behavior  of  optical  systems.  Consequently, 
we  first  devote  several  sections  of  the  paper  to  a detailed  mathemat- 
ical discussion  of  the  polarization  behavior. 

DISCUSSION 


Mueller-Stokes  Matrices  for  Optical  Components 


The  basic  properties  of  the  Stokes  polarization  paranjpters 
have  been  described  In  the  classic  texts  by  Chandrasekhar  and 
Born  and  Wolf.  Futher  properties,  as  well  as  their  relation  to 
the  Mueller-Stokes  polarl^a^l^n  matrices,  have  been  treated  by 
Shurcllff  and  other  authors.  * * 


There  are  three  types  of  optical  components  that  can  change 
the  state  of  polarization  of  an  optical  beam.  These  components 
are  1)  a polarizer,  2)  a compensator  (phase  shifter),  and  3)  a 
rotator.  In  this  section  we  present  the  Mueller-Stokes  matrix  for 
each  of  these  components  and  refer  the  reader  to  the  referenced  works 
for  their  derivation. 
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where  0 £ p(x,y)  £ 1.  For  p(x,y)  - 1 we  have  perfect  transmission, 
while  p(x,y)  ■ 0 corresponds  to  no  transmission  (total  absorption). 
For  example,  p >0  for  a Nlcol  prism  or  Polaroid,  so 
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where  'P  Is  the  total  phase  shift  between  the  orthogonal  axes. 
Specifically,  for  a quarter-wave  plate,  4’  * t^/2,  so 
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where  6 Is  the  angle  of  rotation 
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These  expressions  are  valid  only  If  the  axes  of  the  compon- 
ents are  aligned  with  those  of  the  Incident  beam.  If  the  optical 
device  Is  rotated  through  an  angle  @ then  the  transformed  matrix  Is 
given  by 


M(2e)  - Mj^q^(-20)  M Mj^q^(29)  . 

It  Is  possible  to  construct  an  eigenvector  which  provides 
a simultaneous  diagonal  basis  for  the  Mueller-Stokes  matrices  for 
the  polarizer  and  the  phase  shifter,  but  not  the  rotator.  To  show 
this,  the  Stokes  parameters  are  defined  In  terms  of  the  optical  field 
by 

S-EE*+EE*-I  +I 

0 XX  y y XX  yy 

S - E E*  - E E*  - 1 - I 

1 XX  y y XX  yy 

S-  - E E*  + E E*  - I + I 

2 X y y X xy  yx 

S,  - 1(E  E*  - E E*)  - 1 (I  - I ).  ,,, 

3 xyyx  xyyx  (7) 


The  eigenvector  that  allows  Eq  1 and  3 to  be  written 
In  a diagonalized  form  Is 
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and  the  diagonalized  polarization  matrices  for  the  polarizer  and 
the  phase  shifter  are.  respectively. 


P P 
*^x  *^y 


where  the  tilde,  stands  for  diagonalized  matrix  representation. 
The  relation  between  the  diagonalized  forms,  Eq  9 and  10,  and  the 
nondlagonallzed  Muel ler-Stokes  matrices,  Eq  1 and  3,  is 


M - Mp  M Md 


Mp  - (1/2) 
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The  use  of  a diagonalized  matrix  representation  allows  us  to  write 
the  total  system  matrix  for  a system  with  m polarizing  surfaces  and 
n phase-shifting  elements.  This  Is  found  to  he 
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By  direct  matrix  multiplication,  the  total  system  matrix  is  found  to 
be 

M ■ X 


n Pii-  n n p‘i+  n pJi 

i-i  1-1  1-1  i-» 


0 2 


0 2 


The  reader  should  note  carefully  the  form  of  Eq  19* 

The  ABCD  Polarization  Matrix 

Equations  1,  3,  and  19  are  seen  to  have  the  form 

A B 0 0 

B A 0 0 

M - 

0 0 C -D 

0 0 D C / . (2C 

Thus,  for  the  polarizing  matrix,  Eq  1,  we  see  that  A,  B,  C,  and  D of 
Eq  20  are 


A - (1/2)  (p^+  pj) 
C ■ 2 p p 


B - (1/2)  (p2+  pj) 
D - 0, 


wtiile  for  the  compensator,  Eq  3,  we  have 

A - 1,B“0,C  - cos(|i,  and  D - sin'K 
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Henceforth,  Eq  20  will  be  called  the  ABCD  polarization  matrix,  or 
simply  the  ABCD  matrix. 

We  can  show  that  the  following  equation  holds  for  the 
Stokes  parameters  for  a perfectly  polarized  beam: 

2 2 2 2 

S - S,  +S-  + S-  . (23) 

o 1 z J 

Using  Eq  20  and  23,  we  readily  find  that  the  following  relation 
exists  between  the  A,  B,  C,  and  D parameters: 

2 2 2 2 

A - B + C + D . (2A) 


In  practice,  this  relation  is  very  useful  since  it  obviously  allows 
us  to  find  the  fourth  matrix  element  if  the  other  three  are  known. 
Moreover,  the  relation  can  serve  as  a useful  check  on  all  four 
measured  or  calculated  elements. 


Consider  chat  we  now  illuminate  a system  described  by 
the  ABCD  matrix  with  linearly  +45^  polarized  light.  The  Stokes 


so 


for 

the  incident 

light  is  { 1 

. 0, 

with 

Eq  20  yields 

B 

0 

0 

B 

A 

0 

0 

0 

0 

c 

-D 

1 

^3] 

0 

0 

D 

c 

PH 

1 

f 

1 

» 

j 

c' 

^ j 

1 *3/ 

D / 

o/. 


(25) 


(26) 


COLLETT 


n 


Thus  we  have  the  very  useful  result  that  If  an  optical 
system  consists  only  of  polarizers  and/or  phase  shifters  and  If  the 
optics  are  Illuminated  with  linearly  +45  polarized  light,  each  of 
the  measured  Stokes  parameters  Is  equal  to  one  of  the  ABCD  matrix 
elements. 

The  ABCD  matrix  elements  can  also  be  extracted  using  right 
circularly  polarized  light.  In  this  case,  the  Incident  Stokes  vector 
Is  {l,  0,  0,  l)  , so 
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(28) 


We  note  that  D and  C are  Interchanged  when  Eq  28  Is  com- 
pared with  the  result  for  linearly  +45°  polarized  light,  Eq  27. 

The  previous  analysis  assumes  that  the  axes  of  the  trans- 
mitted field  and  the  optical  system  being  Illuminated  are  aligned. 
In  general,  this  Is  not  true.  Nevertheless,  the  A,  B,  C,  D param- 
eters can  still  be  determined.  This  can  be  done  by  transmitting  two 
sequential  orthogonally  polarized  pulses.  The  optimum  choices  are 
either  and  left  circularly  polarized  light  or  linearly  +45 

and  -45  polarized  light.  The  angle  of  rotation  between  the  axes 
can  then  be  determined  and  the  A,  B,  C,  D parameters  found.  In 
practice,  this  operation  Is  always  done  with  the  polarlmeter. 
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Measurement  of  the  Stokes  Polarization  Parameters  and  the  ABCD  Matrix 
Elements 

We  have  seen  that  each  of  the  Stokes  parameters  is  equal  to 
one  of  the  ABCD  matrix  elements  for  polarizing  and/or  phase-shifting 
systems*  In  this  section  we  discuss  the  measurement  of  the  Stokes 
parameters,  or  equivalently,  the  ABCD  matrix  elements. 

The  Stokes  parameters  of  an  optical  beam  can  be  measured 
by  the  appropriate  positioning  of  a wave  plate  followed  by  a polar- 
izer, set  at  an  angle  @ . This  arrangement  Is  shown  below  In  Fig.  1. 


Figure  1.  Measurement  of  the  Stokes  polarization  parameters. 

A straightforward  analysis  shows  chat  the  detected  Inten- 
sity, 1(q,  (J»),  Is 

1(0,  <l>)  - (1/2)  lSQ+SjC08(2e)  + S2sln(2e)cos(}) -S2sln(2e)8ln(t)]  , (29) 

where  S^,  S^,  S^,  and  S^  are  the  detected  Stokes  parameters,  0 Is  the 
angle  between  ^e  horizontal  x-axls  and  the  polarizer  transmission 
axis,  and  Is  the  phase  shift  between  the  orthogonal  axes  of  the 
wave  plate. 
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The  Stokes  parameters  can  be  obtained  by  sequentially 
setting  0 to  either  0,  7t/2,  or  7t/4  and  Inserting  or  removing  a 
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quarter-wave 

intensities 

plate 
will  be 

in  the  optical  train.  The  respective 

detected 

1(0,  0)  - (1/2) [Sq+S^J , 

(30a) 

I(TT/2,  0)  - (1/2)  [Sq-S^]  , 

(30b) 

I(tt/4,  0)  - (1/2)  [SQ+S2]  , and 

(30c) 

I(Tr/4,  Tr/2)  - ( 1/2)  ISQ-S3]  , 

(30d) 

from  which  the  Stokes  parameters  or  the  equivalent  A, 
parameters  are  found  to  be 

B,  C,  D 

Sg  - 1(0,  0)  + I(TT/2,  0) 

(31a) 

S3  - 1(0,  0)  - I(it/2,  0) 

(31b) 

S2  - 2I(7r/4,0)  - 1(0,0)  - I(7r/2,  0) 

(31c) 

S-  - 1(0,  0)  + I(TT/2,  0)  - I(TT/4,  TT/2) 

(31d) 

In  the  measurement  shown  in  Fig.  1,  the  assumption  was  made  that  the 
source  of  radiation  was  continuous.  This  allows  time  for  selective 
setting  of  the  polarizer/wave  plate.  We  can  represent  the  continuous 
beam  measurement  schematically,  as  shown  in  Fig.  2,  below. 


Incident  Beam 


Wave  Plate/Polarizer 
Single  Channel 


Figure  2.  Sequential  measurement  of  the  Stokes 

Thus  an  incident  beam  propagates  through  a single  optical  train 
(referred  to  as  a single  channel),  and  four  sequential  settings 
and  measurements  are  then  made. 
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It  Is  clear  that  the  CW  method  for  determining  the  Stokes 
parameters  for  a nanosecond  pulse  Is  Inappropriate,  and,  therefore,  a 
different  arrangement  Is  required.  This  can  be  accomplished  by 
Introducing  four  separate  channels,  which  correspond  to  each  of  the 
equations  In  Eq  30  . This  four-channel  polarlmeter,  shown  In  Fig.  3, 
permits  the  simultaneous  measurement  of  the  Stokes  parameters.  A,  B, 
C,  and  D elements. 


Figure  3.  Parallel  measurement  of  the  Stokes 
polarization  parameters. 

A quarter-wave  plate  Is  placed  over  the  fourth  channel,  I(Tf/4,  'n'/2). 

The  arrangement  shown  in  Fig.  3 can  be  schematically 
represented  by 


Incident  Beam 


Figure  4.  Configuration  of  the  four-channel  analyzer. 

We  emphasize  that  this  arrangement  of  the  polarlmeter  has  no  moving 
parts  and,  therefore,  allows  either  continuous  or  pulsed  polarized 
radiation  to  be  measured.  Because  we  actually  use  the  polarlmeter  to 
measure  the  ABCD  parameters,  the  polarlmeter  Is  known  as  either  the 
four-channel  polarlmeter  or  the  ABCD  polarlmeter. 


COLLETT 


In  the  following  and  final  section  we  describe  the  imple- 
mentation of  the  polarlmeter  and  Its  Integration  into  an  automatic 
measurment  system.  This  is  done  using  an  HP  9825A  desk  calculator 
and  a Hewlett-Packard  Interface  Bus  System,  In  conjunction  with 
the  HP  6904B  multlprograouner. 

t 

Computer  Interfacing  in  the  Four-Channel  Polarlmeter 


During  the  past  decade  there  has  been  remarkable  progress 
In  the  development  of  minicomputers,  desk  calculators,  and  Interfac- 
ing equipment.  In  view  of  this  fact  and  the  large  amount  of  data  to 
be  collected  and  analyzed,  the  Hewlett-Packard  982SA  desk  calculator 
was  selected  to  operate  the  four-channel  polarlmeter. 


The  interfacing  between  the  calculator  and  the  polarlmeter 
was  accomplished  using  the  Hewlett-Packard  Interface  Bus  System 
(HP-IB) . Briefly,  the  HP-IB  transfers  data  and  commands  In  parallel 
between  the  components  of  an  instrumentation  system  and  the  HP 
9825A  calculator.  There  are  16  signal  lines,  of  which  eight  are 
data  Input/output  (I/O)  lines  reserved  for  the  transfer  of  data 
and  other  messages  In  a byte-serial,  bit-parallel  manner.  The 
remaining  eight  signal  lines  are  used  for  data  byte  transfer  control 
(3)  and  general  Interface  management  (5). 

In  order  to  provide  even  greater  flexibility,  Hewlett- 
Packard  has  developed  a general  laboratory  Instrument  known  as  the  HP 
69A0B  multiprogrammer.  The  system  Is  placed  In  parallel  on  the 
HP-IB.  This  device  allows  a wide  choice  of  measurements  or  instru- 
ment controls  to  be  made.  The  HP  6940B  multiprogrammer  holds  up  to 
15  plug-in  cards  for  various  purposes;  e.g.,  control  of  stepper 
motors,  voltage  monitoring,  relay  control,  etc.  In  our  system,  an  HP 
693330A  relay  card  was  Inserted  into  the  multiprogrammer  and  used 
to  sequentially  open  and  close  relays  connecting  the  detector/ 
amplifiers  to  a digital  voltmeter. 

While  a voltage  monitor  card  Is  available  for  use  In 
the  multiprogrammer  and  permits  readings  of  150  times  per  second, 
it  Is  limited  to  12-blt  analog  to  digital  conversion.  Since  reso- 
lution Is  of  greater  Interest  than  speed,  an  HP  3490A  digital 
voltmeter,  which  has  a resolution  of  1 microvolt  on  a 100-mllllvolt 
scale,  was  used  on  the  HP-IB.  Because  of  this  much  greater  resolu- 
tion, the  number  of  readings  per  second,  however,  is  reduced  to  5 per 
second. 
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16-bit  format.  In  order  to  make  this  conversion,  an  HP  595000A  inter- 
face unit  must  be  used  between  the  calculator  and  the  multiprogramroer. 

We  pointed  out  earlier  that  the  configuration  of  the 
four-channel  polarimeter  can  be  used  with  a CW  optical  source  as  well 
as  a pulsed  optical  source.  A description  of  the  CW  mode  will  be 
given  first. 

To  Illustrate  the  basic  technique  of  determining  the 
Stokes  parameters,  a block  diagram  of  the  four-channel  polarimeter 
and  the  computer/interf aclng  is  shown  in  Fig.  5.  Several  versions  of 
the  four-channel  optical  head  were  built,  using  Polaroids  initially 
and  precisely  mounted  calcite  prisms  and  quarter-wave  plates  later. 
Each  Polaroid  was  cut  from  the  same  sheet,  and  each  channel  nor- 
malized to  the  first  channel.  The  same  was  done  to  the  detector/ 
amplifier  channels.  The  alignment  of  the  Polaroids  was  measured  to 
be  within  10'  of  arc,  while  the  calcite  prisms  polarizers  were 
measured  to  be  within  2'  of  arc. 


Figure  5.  Block  diagram  of  the  computer-controlled  polarimeter. 

The  four  channels  of  the  optical  head  wer«  uniformly 
illuminated,  detected,  and  converted  from  a current  a voltage 
signal,  using  UDT  lOlA  trans impedance  amplifiers.  B>  programming 
the  calculator,  each  relay  was  closed  for  M seconds,  and  a maximum  of 
five  readings  per  second  was  made  on  each  channel.  This  was  done  for 
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all  four  channels  and  the  measured  data  were  then  stored  in  the 
calculator.  At  the  end  of  the  read  cycle,  the  average  intensity 
reading  In  each  channel  was  determined  and  converted  to  the  Stokes 
parameters,  or  equivalently,  to  the  A,  B,  C,  and  D parameters.  The 
results  were  then  printed  out  on  an  HP  9871A  Impact  printer.  In  some 
Instances  it  was  more  useful  to  display  the  corresponding  polariza- 
tion ellipse,  using  an  HP  9862A  plotter.  At  a future  date  a CRT  will 
be  used  to  speed  up  the  data  presentation. 

The  system  configuration  for  a pulsed  source  Is  similar  to 
that  of  Che  CW  operation.  The  difference  is  the  Increase  In  complex- 
ity of  the  detection  circuit.  Fig.  5,  and  the  addition  of  an  HP 
69434A  Event  Sense  Card  and  an  HP  69600B  Programmable  Timer  Card  In 
Che  multiprogrammer. 

The  detection  circuit  converts  the  nanosecond  pulse 
current  to  a voltage  chat  Is  held  by  the  peak  detector.  The  laser  Is 
fired  once  per  second,  and  a trigger  Is  sent  to  the  Event  Sense  Card. 
The  output  of  this  card  Initiates  the  calculator  for  reading  each  of 
the  relays.  Each  relay  is  programmed  to  be  closed  for  100  millisec- 
onds and  sequentially  read  at  Che  beginning  of  each  200-mllllsecond 
Interval.  All  four  relays  have  then  been  read  at  the  end  of  800 
milliseconds.  During  the  last  200  milliseconds  of  Che  1-second 
cycle,  a single  1-mllllsecond  pulse  Is  sent  to  field  effect  transis- 
tors from  the  HP  6900B  Programmable  Timer  Card.  This  pulse  dis- 
charges the  capacitors  In  the  detection  circuit.  Finally,  after  the 
laser  has  fired  N times,  usually  taken  at  N-20,  the  data  are  averaged 
and  the  A,  B,  C,  D parameters  printed  out. 

CONCLUSIONS 

■ 'V. 

'\ 

We  have  shown  that  a four-channel  polarlmeter  can  be  used 
CO  determine  Che  polarization  characteristics  of  optical  components, 
systems,  and  other  related  phenomena.  While  the  primary  purpose  of 
the  four-channel  configuration  Is  to  measure  the  four  Stokes  polar- 
ization parameters  In  a nanosecond  duration,  there  Is  an  additional 
benefit  to  the  design.  This  Is  shown  by  Illuminating  an  optical 
system  with  right  circularly  polarized  light  or  linearly  +43 
polarized  light.  Each  of  the  measured  Stokes  parameters  Is  then 
equal  to  one  of  the  unknown  A,  B,  C,  D matrix  elements.  In  many 
situations  of  military  Interest  this  property  could  be  Invaluable. 

The  development  of  the  polarlmeter  described  In  this 
paper  evolved  from  Che  recognition  that  all  optical  beams  carry 
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polarization  information.  If  this  information  is  processed  in 
real  time  in  the  manner  described,  then  it  may  be  possible  to  pro- 
vide significantly  Improved  performance  of  military  systems. 
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